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^ . Abstract 

o : 
o 

We complete the formulation of the equations of motion of a non-Abelian 
| gauge field coupled to fermions on a finite-element lattice in four space-time 

c- 1 

dimensions. This is accomplished by a straightforward iterative approach, 

CL' 

qj ' in which successive interaction terms are added to the Dirac and Yang-Mills 

equations of motion, and to the field strength, in order to preserve lattice 

•i-H . 

gauge invariance exactly, yielding a series in powers of ghA. Here g is the 



C3 



coupling constant, h is the lattice spacing, and A is the gauge potential. 
Gauge transformations of the potentials are determined simultaneously. The 
interaction terms in the equations of motion are nonlocal, and can be ex- 
pressed either by an iterative formula or by a difference equation. On the 
other hand, the field strength is locally constructed from the potentials in 
terms of a path-ordered product of exponentials. 
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I. INTRODUCTION 



An intriguing approach to quantum field theories, called the finite element method, has 
been under development for the past decade. (For a recent review see M.) This technique 
of putting the operator equations of motion on a lattice has a number of virtues: 

• It is exactly unitary, or canonical, in the sense that the equal-time commutation rela- 
tions are preserved at each lattice site. 

• It is the most accurate possible method of discretizing the continuum equations of 
motion: Typically, after N steps through the lattice, the relative error is only 0(1/N 2 ). 

• Because it is an analytic technique, it is not subject to the statistical errors of Monte 
Carlo methods. 

• It appears to resolve the fermion-doubling problem, while breaking chiral invariance 
with the expected axial- vector anomaly 0|§. 

• Although a lattice Lagrangian does not exist in Minkowski space, it is possible to 
define a Hamiltonian in the sense of a time-evolution operator, which can be used to 
extract spectral and other information about the theory [£|||. 

Very early in this development it was discovered how to formulate an Abelian gauge 
theory in the finite-element scheme^ @. This was a nontrivial accomplishment, because, a 
priori, it is not obvious that the variational gauge-covariance equations are integrable. It has 
seemed to be much more difficult to formulate a non- Abelian theory. The essential develop- 
ment took place in the late 1980's ||. There it was shown, in two space-time dimensions, 
how to write down the interacting Dirac equation and the Yang-Mills equations, where the 
interactions were determined simultaneously with the gauge transformations. The integra- 
bility of this process was again highly nontrivial. The only thing not completely determined 



1 A rather different approach to the U(l) gauge theory on a finite element lattice is given in jj]]. 
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at that time was the form of the construction of the field strength in terms of potentials. 
(This was given only out to fourth order in potentials.) 

The generalization to four dimensions seemed straightforward. The form for the equa- 
tions of motion was immediately generalizable, but only the first two terms in the field 
strength were given in 1990 ||. Attention shifted to lattice QED 

I returned to this problem in January of this year. I computed the first three terms in the 
field strength in four dimensions, and then recognized the pattern in two and four dimensions. 
The construction is remarkably simple, and should have been recognized very early on: The 
field strength is simply given in terms of a path-ordered product of exponentials of the 
natural finite-element gauge-potential variable (link operators), with further link operators 
providing the necessary gauge transformations in the transverse directions. From the result, 
verification of gauge covariance is immediate. 

The organization of this paper is very simple. In the next section we write down the free 
lattice equations of motion and explain how they are gauged. The full form of the Dirac 
equation, Yang-Mills equations, and the gauge transformation will be given. In Section III 
the results of the iterative calculation for the field strength through order g 2 hA 3 are given, 
where g is the gauge coupling constant, h is the lattice spacing, and A is the vector potential. 
In Section IV we infer the exact form of the field strength, which is immediately verified 
to be gauge covariant. The continuing direction of this research program is sketched in the 
concluding Section V. The essential element of the proof of gauge covariance is given in the 
Appendix. 

II. FINITE-ELEMENT EQUATIONS OF MOTION AND GAUGE INVARIANCE 

The two-dimensional work has been described in detail |]T]||, so we will begin immediately 
in four- dimensional Minkowski space. Let us start from the free Dirac equation in the 
continuum, 

^ 7 ^ + /i)V = 0. (2.1) 
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The linear finite-element prescription consists in replacing derivatives by forward differences, 
while in directions in which there is no differentiation, a forward average is employed. If we 
use the following notation for forward averaging, 

•Em ~^{%m+l ~\~ •^rn)i (2'2) 

the free finite-element lattice Dirac equation is 

-j-{iha,n+i - iha,n) + -^■(^mj+^m^n ~ i> mj ,Ta ± ,n) + ^m,n = 0, (2.3) 

where h is the temporal lattice spacing, A is the spatial lattice spacing, m represents a 
spatial lattice coordinate, n a temporal coordinate, and the _L indicates directions other 
than the one singled out. 

Gauge transformations are introduced just as in the continuum. The Dirac equation 
( p.3| ) is invariant under the infinitesimal global phase transformation 

#0m,n = ig5u4> m , n - (2.4) 

However, if Sou depends on position, invariance is spoiled. Let us define a local transformation 
so that the mass term in ( |2.3p transforms covariantly: 

8ipm,n = igSoj m , n ihi,n- (2.5) 

The forward differences in ( |2.3|) do not transform covariantly; to achieve covariance we must 
add interaction terms. The interaction is with a vector potential transforming by the lattice 
analog of 

A^A^ + S^A^, 6^A^ = d^, (2.6) 

that is, 

5 (0) (A)) m = ^(5Anr, n+ i - <5A m , n ), S^iA^ = -r(5A mj+ i >m± ,n - SA mj ,m±,ri), ( 2 -7) 

where m represents the four- vector coordinate (m, n) and the superscript (0) is a reminder 
that this is the first term in a series of variations. We will take the connection between 5u 
and SA to be the finite-element one: 



Su m = 5A m . 



(2.8) 



Then, using the boundary conditions (where M is the number of lattice sites in any spatial 
direction) 



VW,-+M,m_L,n — ( 1) + 1prrij,m ± ,n-, $' l Prnj+M,mj_,n — ( 1) + $1prnj ,m ± ,ni 



(2.9) 



and the initial condition 



84>m,0 = y(^m,0 + Su m -OV^n.O, 

it is easy to see that the noncovariance of 

vy , vy* 

7 (^m,n+l V'm,™) "I - ^ (iprrij +I,i5x,n ^Prrij ,m±,n) 

will be cancelled if the following interaction terms are added to the left side of (|2.3|) : 

m-j+M 



(2.10) 



(2.11) 



4 1} = ^A £ (-l) mj+m H^)^, m± >m ; , ni± , n 



(2.12) 



and 



- o n 

K$> = -2^^^'(-l) n+n '(^o) m ^ 



(2-13) 



n'=0 



where the prime on the summation sign indicates that the n' — term is counted with half 



weight. The variations to which the interaction terms are subjected are, from ( p. 71) 

^ ^ iA-j)mj— l,m_L,n ~~ T(^mj,mi,)i ~~ &^mj— l,mj_,Ti) 

^ (0) (^o) m ^ = T(^m,n - <^m,n-l)- 
ft 



(2.14a) 
(2.14b) 



It will be noted that the form of the interaction term with the scalar potential involves only 
the fields at the current and earlier times. This feature persists, and enables one to solve 
the equations by time-stepping through the lattice. It is easy to verify |6]|| that the Dirac 
equation with fl2.12[ ) and ( |2.13p included reduces to the correct continuum limit. [In this 
verification, the form of the initial condition ( |2.10|) is crucial.] 



However, we are not finished, because we have not varied the interaction terms (|2.12|) 



and ( |2.13| ) with respect to ip, ( |2.5| ), nor achieved covariance of the interaction term. This 
can be accomplished by adding new interaction terms 1^ and K^ 2 \ and additional gauge 
transformations, denoted by A straightforward calculation reveals that 

V<? = to***® - 8 W l£> - 5 (0) 4 2) . (2.15) 

Here the first term on the right is the required covariance term, the second term involves a 
new variation of A, 

and the third term is the S^ ' variation of a new interaction term 

m' =m 7 +l m'!=m'-+l 



Similarly, we find the new variation of A 



0) 

n>l: 5 (1) (^o) m ,^r = -w [^m,n + ^m,n-i, (A)) m — ], (2.18a) 
n = : 5 (1) (^o) m ,3T = y[^w m , + &*> m ,-i, (A ) m ^— } 

+ j[Su mfl - 6u m -!, (A ) m — ], (2.18b) 

and a new interaction term involving the scalar potential, 

. q n n' 

K$ = -^{ighf Y:'^ l ) n+n '^m— (A ) m —^m,n». (2.19) 

' n '=0 n"=0 

Here the double prime on the second summation sign means that both the first and last 
terms are counted with half weight, and that if the upper limit is zero, the sum vanishes. 
Of course, the variations (|2.16 ) and ( |2.18[ ) precisely reduce to the usual non-Abelian 



transformations of a gauge field in the continuum, and the new interaction terms ( |2.17| ) and 
( |2.19| ) vanish in the continuum limit. 

This iterative procedure continues indefinitely. Fortunately, it exponentiates. Perhaps 
the easiest way to express I^ N \ where N is the order in gAA, is by an inductive formula. 
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Let us adopt a shorthand notation for ease of writing the results. Let I = and because 
Ij depends locally on all coordinates except rrij, let us suppress those other coordinates, call 
rrij = m, and simply write I m to stand for (Ij) mj , m± ,n- Then we find 

1 m+M N 1 

N>1: I^ = l - (-l^E^Mrf)^, (2-20) 

m'=m+l fc=l 

where 

"4m = (^4j)mj-l,mj_,n (2-21) 

and we define = — 2i-0 m .. i5ix ,n/A. The gauge transformations of the potentials are given 
by nested commutators: 

k^l: 8^ Am = ^^-yf [• ■ • [^m - <W-i, AJ, • • -,A m ], (2.22a) 

fc nested commutators 

S (1] A m = j[5uJ m + 5u m -l,A m ], (2.22b) 
where Ek is the fcth Bernoulli number. The required covariance statement 

N 

V£° + E ^ (fc) e- fc+1) = igfami™ (2.23) 

fc=0 



is proved in the Appendix of || . By summing ( 2.20| ) over all N we derive the sum equation 



satisfied by the full interaction term involving the vector potential, / = ^at=i I ■ 



m+M 
m'=m+l 

from which a difference equation can be immediately derived: 

AA 2^ ■ a.4 

Im ~\~ & m Im— 1 = ^ (l 6 3 m ) Ipmj ,m ± ,n- (2.25) 

Of course, I transforms in the required manner (5 a = Y2T=o 

5^I m + 5 A I m - 5 (0) 4 1} = igSu m I m (2.26) 
where the third variation on the left-hand side is just that necessary to cancel the nonco- 



variance of the free Dirac equation ( |2.3|) . 



In just the same way we can derive analogous expressions for the full interactions and 
transformations of the scalar potential. Again we suppress all but the time coordinate n, 
and write the interaction term as K — > 7°A" in what follows. The only complication has to 
do with the initial conditions. The order N interaction term is given by 



N 



n'=l 



k=l 



-Wh) k k {N -k) {igh) N 



-1)VC^ 0) , 



(2.27) 



where Kn^ = —2itp m:n /h. The required transformation law 

JV 



(2.28) 



fc=0 



is satisfied provided the transformations are given by 



k^l: 5 (k) A r , 



h k\L 



i "4-n]> 5 «^rj) ^ 7^ 0, 



k nested commutators 



5^ An = l jr[5uj n + 6u n -i, An], n^O, 



k^l: 5 {k) A Q 



-igh) k B k 
h 2 k k\ 



[■ ■ ■ [5uj - 5u_x,A ], ■ ■ ■ , Aq], 



k nested commutators 



5 {1) A = y[^o + 6u}-.x,Aa] - ig^idwo - Su^Aq] 



(2.29a) 

(2.29b) 
(2.29c) 

(2.29d) 



A sum and difference equation may be immediately derived for the full interaction term 
involving the scalar potential K n = Y1n=i Kn '■ 



K n = - J2(-l) n+n ' {e- t9hA "> - 1) (K n , + + (-1)" (e l9hA °/ 2 - l) , (2.30) 



n'=l 



and 



1i 

K n + e^Kn-x = - (1 - e i9hA -) fe, r 



(2.31) 



The full lattice Dirac equation is given by (|2.3j ), fl2.20|) , and ( |2. 271 ): 
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h 



(^m,n+l - lfei,n) + K v 



A 



(VWj+l,m_L,n VW,- ,mj_ ,n) \Ij 



+ /^0m,n = 0. 



(2.32) 
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We have already emphasized that this equation gives ip m ,n+i in terms of fields at time n and 
earlier, so that this difference equation may be solved by time stepping through the lattice. 

We have given so much detail on the Dirac equation because this is the hard part of the 
calculation. The crucial point is that we have now determined not merely the interaction 
terms in the Dirac equation, but the precise form of the gauge transformations on the lattice. 
We must now use that information to construct the gauge sector of the theory. (If one starts 
directly with the gauge sector, it is much less clear how to proceed.) 

Let us assume that we can construct a gauge-covariant field strength from the poten- 
tials. (That construction will be given in the next two sections.) That is, we have a field 
strength defined on the lattice which transforms under an infinitesimal gauge transformation 
according to 

S{F^) m = ig[5u m , (F^) m }. (2.33) 



(This is precisely analogous to ( |2.5| ).) This is a difference equation for the transformations 
of the field strength at the lattice sites. Because we require boson fields to be periodic, we 
must henceforward assume that M is odd. (Recall the periodicity condition (|2.9j).) For the 
inital condition we adopt the analog of ( |2.1(J| ), 



<5(-Foi)m,o = — [&L>m,o + Su m - h (F i)m,o\- (2.34) 

Our starting point for describing gauge fields coupled to a fermionic current is the "free" 
Yang-Mills equation 

T } 1 [(-^ )mu+l,m f i,m± ~ {F^ ) m„ ,m M ,m jj = (j'^)m- (2.35) 

Here, for ease of notation, we have not destinguished between the space and time lattice: h 
or A should be used as appropriate. A gauge covariant candidate for the fermionic current 
is (T is the generator of the gauge group) 

# = ^ryvw (2.36) 
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because, then, according to (|2.5|) , the right side of ( |2.35| ) transforms covariantly: 



6& = ig[6u m ,j£\. (2.37) 

It is the difference in ( |2.35| ) that spoils gauge covariance. But now it is immediately obvious 
that the previous analysis for the Dirac equation applies. Multiplication simply gets replaced 
by commutation. We simply must add to the left side of ( |2.35|) the interaction term 

JM = ^/c-, /C£ = /C (2.38) 

(that is, in K? the only nonlocal dependence is in m v = n), where JC U satisfies the difference 
equation 

JC n e^ hA - + e^/C^! = ~ [e i9hA \ (F^) m „ >ffij .] . (2.39) 

The solution to this difference equation can be immediately written down by comparison 
with the corresponding terms in the Dirac equation (multiplication being replaced by com- 
mutation) and is given explicitly in [IjlEJ. In particular, the inductive formulae are 

i rn+M N , . .s fe 

ICgP = \ E (~ 1 ) m+W E L 4r L J^- ^ L_J ^ ^" fc) ] • • ■]]]; (2-40) 

m'=m+l k=l , , 

fe nested commutators 

for the v — j interaction, and 

N , . Mfc 



n'=l fe=l ' ' S ~^ 



(N-kh 

— I . /-\ . . I . I . /-\ . . I . I * * * . I . r-\ . . I . !\ , 

k\ 



k nested commutators 

+ WM^ l^l [ ±1 •••]]]; (2-41) 

N nested commutators 



for the v — interaction. 



III. ITERATIVE CONSTRUCTION OF THE FIELD STRENGTH TENSOR 



The finite-element prescription for the construction of the field strength in terms of 
potentials, in lowest order, is simply the lattice curl: 
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{FpJ)m — , [(^)m M +l,m ± (^4.i/)m M ,mj_] ■ 

(3.1) 

(In this section, for simplicity, we set h = A; the distinction can be readily recovered if 
desired.) Of course, this is invariant under the transformation S^°\ ( f2.14j ). Next, we apply 
the transformation (|2.16|) and (|2.18|) ; remarkably enough, we are able to write the result 



in the form 

^ViJV = ig[5u m , (F<$U - 5(°\F$) m , (3.2) 



that is, the difference from the required covariance fl2.33j ) is a total variation of a new 
interaction term. To present this interaction most simply, let us introduce some further 
notation. We use the averaged potential (|2.21| ), which we write here as 



•A-klmn ~ 2 i^klm-a + ^klmn-l) > (3-3) 

where henceforward we order the coordinate indices k, A, fi, v. (Here n = m v is not 
necessarily the time coordinate.) In terms of this, we define a "field strength" variable local 
in the transverse coordinates (k, A) 

klmn , [(A (A vjklmn+1 

(A) 

klm+ln+1 klm+ln 

] . (3.4) 

Under the <5 (0) variation (|2.14j) , 

6^°\A u )klmn = -ri^klmn ~ ^klmn-X)i (3.5) 

h 

Tp) is invariant. 

In terms of these new variables, the lowest-order field strength is obtained simply by 
averaging in the transverse coordinates: 

klmn' 

(3.6) 



Because the first-order variation is given by fl2.16|) and ( |2.18| ) (except for Aq at the initial 
time), 
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klmn klmn 

I (3-7) 



the interaction term found in (13.21) can now be written as 



i^jlu )m (3~jiv ) klmn ^ {[(A K )k+llmn, (^)k +llmn ] + [(Ax) kl+lmn, fci+lmnj 
+ 2 [(Ak) k+ll+lmn, (^)fc+ll+lmn] + g [(40fc+lH-lmn> (^"j^}) k+ll+lmn] 



+ 2 [(A) (^)wWmn]}) (3-8) 



where 



klmn „ -{[(^ fijklm+ln, 

(4) 

klm+ln+1 ]-[(A v jklmn+1 i (A v)klm+ln+l\ 

- [(4) 

klm+ln+1 1 (A v ) klm+ln+1 ] + [(A fijklm+ln, {A v ) klmn+1] 

- [(4) 

klm+lni 

{A 

v ) klm+ln+1 ] - [(4) 

klm+ln+1 1 (A u ) Mmn+l]} ■ (3.9) 



It will be noted that ( |3.9| ) is exactly the form of the first-order field strength found in two 
dimensions ||. 

We continue, by finding the second-order field strength from the covariance equation 

£ (2) (^k + 8<»(F<2) m + 5^(F$) m = ig[8u w (F^U (3.10) 

Here (n > 0) 

g 2 h 

v ) klmn --T^-[[8Uklmn — SLUklmn-li (A v ) klmn \,(A v )klmn\, (3.11) 



from ( 2.22|) and (|2.29 ). The result of a long calculation is remarkably simple: 



(F^J)m {J~ l J) klmn ^ {[(A k) k+llmni k+llmn kl+lmni (^) kl+lmn\ 

+ ~[(A k) k+ll+lmn, (^)k +ll+lmn } + t;[(Ax) k+ll+lmni (^J)k +ll+lmn\ 
Klk+llmn, (^) k+ll+lmn } + t;[(Ax) kl+lmn, (^u) k+ll+lmn]} 



kl+lmni [(Ax) kl+lmni kl+lmn\ 



[(A K ) k+llmn i [(Ak) k+llmn: (^jiJ)k+llmn\\ 
+ 2 [(A K ) k+ll+lmn, [(A K ) k+ll+lmn, (J~$)k+ll+lr 
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k+ll+lmni k+ll+lmni {^~ji,J)k+ll+lmn]\ 

Kik+llmnt [(A Kjk+llmni (^) k+Yl+Ymn\\ 



kl+lmni Ux) kl+lmni 



kl+lmni ll+lmn i 



) k+ll+lmn\ 



\\Ak) k+llmni 



where 



(F$hmn = -§{[^i> [Bi, C x \] - [Cx, [&, BJ] - [B u [B x , Co]] + [d, [d, B ]] 
+ 2[B 1; [B , d]] - 2[d, [Co, Bi]] + 2[B , [B u d]] - 2[C , [d, Bi]] 

- 2[B 1; [B , Co]] + 2[C 1; [Co, B ]] + 4[B , [B l5 C ]] - 4[C , [d, B ]] 
+ [B , [B , d]] - [Co, [Co, B x \] - [B , [B , C ]] + [C , [C , B ]] 

- [B , [B , B x }\ + [Co, [Co, d]] - [Bx, [B , B 1 ]] + [d, [C , d]]}. (3.13) 

Here we have adopted one final bit of condensed notation: 
B x = {AJ klm+ln+1 i B = (A) 

klm+lm v ) klm+ln+l i v ) klmn+1 ■ 

(3.14) 

One additional explicit term in JF, for J/w , was given in [§]. 

It is quite clear that this process may be continued indefinitely. We will now see that 
these terms may be summed exactly. 



IV. EXACT FIELD STRENGTH 

At first glance, the pattern for the construction of the = Y1n=i seems a bit 
mysterious. However, upon a bit of inspection, the terms given here and in || are entirely 
reproduced by the expansion in h of the following path-ordered product of exponentials (link 
variables) around the finite element (plaquette): 

e -ift 2 g(^„) fc!mn _ p e ~igfA-dl _ e ~ihgB e ~ihgCi e ihgB! e ihgC (a ]\ 
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This is precisely the form conventionally used in lattice gauge theory (implicitly in Wilson's 
original paper |IT| , and explicitly, for example, in fll2|). Of course, here, the appropriately 



averaged finite-element potentials are employed.^ That this is exactly correct emerges im- 
mediately when we recognize the effect of the gauge transformations on the link operators 

5e -igh(Au) k i mn+ i = ig {5Lu klmn e- i9h{A » )klmn+1 - e- i9h{A » )klmn+x 8u klmn+1 } , (4.2) 

and the corresponding adjoint equation. The proof is given in the Appendix. It then follows 
immediately from ([4.1|) that 

}. (4.3) 

Now, we immediately recognize the structure of the full field strength = X]jv=o Fffi 1 
the first three terms of which are given in fl3.6| ), ( |3.8| ), and ( |3.12| ): 



kl-\-lmn 

y 1 HvJklrmi ~ / ^^\- r nv jklmn T e \~r fiu J kl+lmn^ 

_|_ e -igh(A K )k+iimn (^p ) fc+i; pigHA^k+iimn 



1 p-igh{-A\)kl+lmnp-igh(A K )k+ll + lmn ( -T \, , pigh(AK.)k+u+lmn p igh{A\)kl+Vm.n 

~T~ 2 c c \ J l-LV ) fe + lt+lmn c c 

limn p -igh(A x ) Jgh(Ax) igh(A K ) k +u run 

\ 2^ c \V fxv ) ft+li+lmn c c 

(4.4) 

That i 7 ^ transforms covariantly according to ( |2.33| ) follows immediately from ( |4.2[ ). That 
is, the field strength is the average of the "transversely-local" field strength over the 
four corners of the transverse finite element, with gauge factors inserted appropriately to 
gauge-transform the field strength back to the origin (klmn). 

Finally, we note that at the initial time F uv involves (^4o)m,i and not (^4o)m,o, so that 



the initial variations Q2.29cj) and (|2.29d| ) never come into play. 



2 That this is nontrivial can be seen from the remarks in [12]: "This identification of A a (n) with 
the vector potential is only strictly correct in the continuum limit," and that only as the lattice 
spacing goes to zero do the link variables transform equivalently to 
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V. CONCLUSIONS 



We have thus finished the task of giving a complete set of gauge-invariant equations of 
motion for Yang-Mills fields coupled to fermions. This should enable us to apply finite- 
element methods of solution which have so far been applied to scalar field theories ]n|] and 
to electrodynamics PH|2]J3fl. New methods are being developed as well ||. It is therefore 
reasonable to hope that we may be able to use this technique to shed new light on difficult 
problems in non-Abelian theories such as QCD. 
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APPENDIX: PROOF OF GAUGE TRANSFORMATIONS OF LINK OPERATORS 



The proof of (|4.2|) follows essentially from the proof of covariance given in the Appendix 
of M . We have, for example, from (|2.22| ) (we are using the notation ( |3.14| ) and only display 
the values of the /i, v indices) (n > 0) 

5 (k) e i 9h B = j2 Mf. . . [6uJm+ln _ 5 Bo ],..., b ] b»- 1 -\ k * 1, 

N=0 ' 1=0 ' N I nested commutators ' 

(Al) 

and 

6 (i) e i 9h B = J- B l l j[5uj m+ln + 5u mn , B }Bt l -\ (A2) 

N=0 ' 1=0 

so the total variation of the link variable is 
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Se ighB = J2 s (k) e i 9 hB = ig [ SuJrnn ,e i9hBo } 



k=0 



oo N-l 



+ T E E B o { -^-^l- • • - ggk^ gol^"'" 1 (A3) 



ft 



AH fc! -t 



fc=0 AT=1 Z=0 

The nested commutator is simply given by 



k nested commutators 



[■■■[6u,B),...,B) = J2(- 1 ' 

v .. / z — ' 



k nested commutators 



s=0 



B s 5uB k ~ s . 



Then the sum on the right hand side of ( [A3j ) reduces to 

j oo r 

- 22 z2(wh) r+1 B%(8uj m+ i n - 5u mn )B t 



yr-p 
5 



r=0 p=0 



X 



E-d 



fc=0 



fc! (r + 1 - fc)! 



E 



s=p— r+k 



V) 



Se l9hBo = ig(6u m+ln e i9hBo - e i9hB °6u mn ), 



(A4) 



(^ m+ln -^ m „)e l9/lBo , (A5) 



where the last two summations over s and k are carried out using identity (A3) of H. Thus 
we have established the desired result: 



(A6) 



which is the adjoint of (|4.2|) . 
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